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In this work we focus on the dynamics of macromolecular networks formed by end-linking identical poly- 
mer stars. The resulting macromolecular network can then be viewed as consisting of spacers which connect 
branching points (the cores of the stars). We succeed in analyzing exactly, in the framework of the generalized 
Gaussian model, the eigenvalue spectrum of such networks. As applications we focus on several topologies, 
such as regular networks and dendrimers; furthermore, we compare the results to those found for regular hy- 
perbranched structures. In so doing, we also consider situations in which the beads of the cores differ from 
the beads of the spacers. The analytical procedure which we use involves an exact real-space renormalization, 
which allows to relate the star-network to a (much simpler) network, in which each star is reduced to its core. It 
turns out that the eigenvalue spectrum of the star-polymer structure consists of two parts: One follows in terms 
of polynomial equations from the relaxation spectrum of the corresponding renormalized structure, while the 
second part involves the motion of the spacer chains themselves. Finally, we show exemplarily the situation for 
copolymeric dendrimers, calculate their spectra, and from them their storage and the loss moduli. 

PACS numbers: 



I. INTRODUCTION 

A series of recent works focuses on the synthesis of macro- 
molecular networks by end-linking star polymers (J Q, 0. EL 
0,0]. Now, previous work |7] has considered hyperbranched 
polymers, which were iteratively built from rather small units: 
there the use of regular repeated units allowed a complete, an- 
alytical determination of the eigenfrequencies. Here we ex- 
tend this work by focusing on structures built from stars. Us- 
ing stars as building blocks allows to vary both /, the func- 
tionality and also k, the numbers of monomers along the arms. 
The stars will be envisaged as being end-linked in such a way 
that two end monomers belonging to two different stars get 
connected by an additional link. In this way we obtain a sys- 
tem of branching points connected by spacers; this system of 
branching points may, on the other hand, be rather arbitrary, 
and may by itself form a complex network such as a regu- 
lar dendrimer Q|, a branched polymer 0] or a branched 
copolymeric structure 1 8 ] . 

The simplest method to determine the dynamics of such 
networks is to treat them as generalized Gaussian structures 
(GGS) l9l fiol fill O : the polymer is viewed as consisting 
of beads that are linked by harmonic springs. In this way, the 
classical Rouse-model gets adapted to general topologies. As 
it is well-known, the Rouse-model describes a chain molecule 
immersed in a viscous solvent, so that for each monomer the 
damping due to the solvent and the harmonic potential due 
to the other monomers are accounted for. On the other hand, 
the Rouse-model neglects hydrodynamic interactions, inertia 
terms, the excluded volume and entanglement effects. Now, 
the limitations of the Rouse model are well understood; al- 
though not very realistic, Rouse-models and their extensions 
are certainly very useful in unraveling the dynamics of com- 
plex structures. I^.fl^l 
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In the study of networks built from stars we will take all 
the springs to be equal, but we will differentiate between the 
monomers at branching points and the monomers along the 
spacers; we will let the two kinds of monomers have differ- 
ent mobilities i.e. different friction coefficients. The reason 
for this are twofold: first, in many cases stars are prepared 
by having a central monomer which differs from the other 
monomers; second, mathematically this generalization does 
not significantly burden the formalism. This procedure, on the 
other hand, allows us to discuss easily the case of copolymers. 
Now, as was shown in former work I13lfl4ll . the relaxation of 
copolymeric structures containing units of different mobility 
often differs qualitatively from the relaxation of similar ho- 
mopolymeric structures, in which all the beads are the same 

id in. 

In what the relaxation dynamics of the networks is con- 
cerned, we focus on the complex dynamical shear modulus, 
G*(u), or in more conventional forms, on its real and imagi- 
nary parts, which are the dynamic storage and the loss mod- 
ulus, G'{u>) and G"(w), respectively. The advantage of look- 
ing at these quantities is, first, that they are widely measured, 
and, second, that there are intimately related to other dynamic 
quantities, such as the dielectric and the magnetic relaxation. 
Moreover, in the evaluation of G'{lu) and G"(lo) only the 
eigenvalues of the underlying connectivity matrix enter, but 
not the eigenvectors 19, M (ill 211 . This fact considerably simpli- 
fies the mathematical treatment, by allowing us to focus on the 
eigenvalues. As we proceed to show we can relate the eigen- 
values of our system of end-linked stars to the corresponding 
underlying spacer-free structure. 

The basic realization here is that the problem considered 
here admits an exact real-space renormalization. This may ap- 
pear surprising after the efforts displayed up-to-now 1 151 Hal , 
which lead only to approximate results. The interesting fea- 
ture is that the previous works used models in which some 
beads were assumed to be fixed 11151 1 1 611 - Surprisingly, it is 
just this assumption which complicates the analytical solu- 
tion; by allowing all beads to move freely the problem sim- 
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plifies (of course, some normal modes leave groups of beads 
immobile, but this is a result, and not an a priori assumption). 
The situation here parallels that found in Ref. |7], where reg- 
ular hyperbranched polymers were investigated; the exact re- 
sults in |7] followed earlier works, which were unable to solve 
the problem exactly, again since they took the positions of the 
peripheral beads to be fixed. 

The paper is structured as follows: In the next section we 
present the model considered and recall the forms of G'(oj) 
and of G"(oj) in the GGS-formalism. In Sec. 3 we explain 
our renormalization procedure, which allows us to relate an- 
alytically the spectra of the networks under investigation to 
those of simpler structures, in which the spacers and the dan- 
gling chains are eliminated. In Sec. 4 we display the use of 
our method in several cases and focus on star polymers, on a 
cube polymer and on topologically periodic networks. Den- 
dritic structures made out of copolymeric stars are analyzed 
in Sec. 5, where we present numerical results for G'(uj) and 
G" (uj). Our conclusions are summarized in Sec. 6, while the 
mathematical details are collected in four Appendices. 



II. THE MODEL 

In the following we consider a system of N end-linked 
polymer stars. Each star consists of a core (a central /- 
coordinated bead) from which emerge / chains of k beads 
each. The central bead may differ (chemically and function- 
ally) from the other beads, a fact which we take into account 
by assigning it a different friction coefficient. In Fig. ^ we 
show such a network formed from end-linked stars. Corre- 
sponding to the coloring used in Fig. ^ we call in the follow- 
ing the core bead "black" and the chain beads "white". Note 
that in Fig. [2 the stars are linked at their ends, through the 
means of I additional bonds, to form a single, connected clus- 
ter. The linking topologies which we consider are rather arbi- 
trary; the restriction that we impose is that the cores (branch- 
ing points) form a simple graph structure. By this we forbid 
links between chains belonging to the same star as well as 
multiple links between two stars. 

Now we model the system of Fig. [T]in the framework of 
generalized Gaussian structures (GGS), which are an exten- 
sion of the Rouse-model for Gaussian chains 

G3IIH0. We 

focus on networks of beads connected by harmonic springs 
with elastic constant K, and we take the networks to be im- 
mersed in a viscous solvent, where the inertial forces can be 
neglected. The corresponding Langevin equations of motion 
for the system have the form: 



dt 



[rijt)-^)] =F(r J - ,*)+£,•(<). (1) 



Equation {0 expresses the fact that for each bead, say bead j 
at position r 3 (t), the friction force and the elastic forces acting 
on it exactly balance the external and the fluctuating forces. 



In Eq. 0, the r io (t) with a = 1, . 



denote the positions 



of the rij neighbors of bead j, Q is its friction coefficient, 
and fj (t) and F(rj , t) are the respective random and external 




FIG. 1: Network consisting of end-linked polymer stars, of which 
one is highlighted by a circle. Here N = 10,/ = 4, and k — 2, so 
that Nf u u = 90. Note that the structure is very flexible so that in 
many realizations it will appear considerably less ordered than here. 



forces; the fluctuating forces are taken to be Gaussian, with 
(fjAtffiA?)) = 3-k B TC 3 8 l0 8 ab 5{t-t'). Moreover, we take 
the friction coefficients of the white beads to be Q = (o and 
of the black beads to be Q = a Co- 
in the GGS picture G'(u>), the storage and G"(uS), the loss 
modulus are given in terms of the eigenvalues \ of the homo- 
geneous version of Eq. Q (where the right-hand-side is put 
equal to zero). One has namely, in a reduced form llOt flTII 



and 



G'(Q) = 



G"(Q) = 



N 



full 



N fu „ 

E 



o> 2 + 4A2 



(2) 



full 



2u>X r , 



N }u ii ^ ui 2 + 4A, 



(3) 



In Eqs. (0 and (0 the reduced frequency Q = ui/tq is given 
in terms of the fundamental time tq = Co/ K, Nf u u is the total 
number of beads and the A m are the Nf u u — 1 non vanishing 
eigenvalues of the system, where 



N full =N(fk+l), 



(4) 



given that each star is composed of a core (black bead) and 
of fk white beads. As shown previously, 1 1 3l 1 1 411 the general 
forms of G'(oj) and G"(u>), well-known from the study of 
homopolymers, also hold for heteropolymeric networks, such 
as the ones considered here. 

To determine the eigenvalues of the left-hand-side of Eq. 
Q we focus on its normal modes, which have the form 



,(m) 



(*) 



exp(- 



A m i/r ). Inserting them in the left- 



hand-side of Eq. Q leads to the following system of equa- 
tions: 



c. 



• 3 \ 

S>0 



(m) 



a=l 







with 1 < j < N 



full- 



(5) 
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With Nfuii being the total number of beads in the system, 
this represents Nf u u equations. As is well-known from the 
general theory, the linear system given by Eq. (0 leads to 
Nf u u — 1 positive eigenvalues, A TO > 0, and to a single van- 
ishing eigenvalue, Ai = 0. 

In the following we look separately at the system of equa- 
tions belonging to each m and for ease of notation drop the 
index m. For the /-coordinated core (black bead) of each star, 
Eq. reads: 

/ 

(f-aX)<Z>j = X)*i-- (6) 

a=l 

For the white beads of the spacers (these are doubly coordi- 
nated), the corresponding equation is 

(2- A)$ i+1 =$ 4 + $ l+2 , (7) 

where we numbered the beads consecutively along the spacer. 
Again from Eq. we have for every peripheral bead k 
(which is sm<^?/-coordinated) 

(1-A)*fc = * fc _i. (8) 

Remarkably, the system of Eqs. 0, and can be sim- 
plified analytically for any arrangement like the one in Fig. [2 
This fact extends to stars the exact mathematical treatment de- 
veloped for small building blocks, used in the study of regular 
hyperbranched structures. 

III. RENORMALIZATION PROCEDURE 

The following section is devoted to the simplification of 
the system of Eqs. 0, and 0. The basic idea is that 
one has an analytical procedure at one's disposal, by which 
one can "delete" most of the beads of the original system, 
namely those on the spacers and on the dangling chains; in 
this way one is led to a smaller, reduced system, consisting of 
the (directly-connected) black beads only. The procedure is 
exact and we will develop it in the main body of this section. 
The method parallels that of determining the eigenvalues of 
regular hyperbranched structures, a problem treated in Refs. 
IU5L Ha fl8l fl9ll . but it is more extensive. Moreover, it will 
turn out that the eigenvalues of the original system fall natu- 
rally into two classes: The first class is the one related to non- 
vanishing eigenvalues of the reduced black system, while for 
the second class the corresponding eigenvalue of the reduced 
system is zero. As a consequence, we follow a step-by-step 
development and focus first on the spacers of a system such as 
displayed in Fig. [2 i.e. on the chains of white beads. 

A. General framework 

Consider now Fig. ^ in which for each pair of linked stars 
the core black beads are connected by a chain of s — 2k white 
beads, and in which each dangling chain consists of k white 



beads. We consider first an arbitrary chain of white beads and 
number them sequentially denoting the starting black bead by 
"0". In this way the amplitudes (eigenmode components) of 
the white beads get denoted by $i, $2> • • •> while that of the 
black bead gets denoted by $o- F° r a dangling chain the last 
amplitude is while for a spacer it is $ s , with $ s +i denot- 
ing the amplitude of the next black bead. 

Now, all the white beads fulfill Eq. 0. Introducing the 
notation 

2x = 2 ~ A, (9) 

Eq. reads exemplarily for the white bead i + 1: 

= 2x$ i+1 - $ i+2 . (10) 

The following considerations are significantly simplified by 
making use of the Chebyshev polynomials U n (x) of the sec- 
ond kind, see Eq. (22.2.5) of Ref. 1 20] . The first two Cheby- 
shev polynomials read: Uo(x) = 1 and U\(x) = 2x, Table 
22.5 of Ref. 12011 . where they are standardized by requiring 
that U n (l) = n + 1. Furthermore, the U n (x) obey the recur- 
sion relation: 

U n {x) = 2xU n - 1 {x) - U n ^ 2 (x), (11) 

see Eq. (22.7.5) of Ref. lEoll . In terms of trigonometric func- 
tions one also has Eq. (22.3. 16) of Ref. 1 20]: 

U n {cosip) = ; , (12) 

smip 

where for (p £ R one has — 1 < x < 1. Of course, as polyno- 
mials the Uf tS (x) are defined over R, which means that (p in 
Eq. dl2l is allowed to be complex. 

We now proceed to show by induction that the amplitudes 
of each pair (n,n + 1) of white beads along the chain are 
related to that of the black bead through 

$o = U n $ n -U n - 1 $ n+l , (13) 

where here and in the following U n (x) = U n is implied, if 
not indicated otherwise. 

Now, for n = l Eq. Jl 3i reproduces Eq. dlOt for i — 0. Let 
us therefore prove the induction step. Inserting Eq. JlOi into 
Eq. dl 31 for n = i we find 

$o = U n (2x$ n+ i - $ n+2 ) - ?7„_i$ n+ i = 
= (2xU n - U n -i)$ n+ i - UA+2 
= U n+ i$ n+1 - t/A+2, (14) 

where the recursion relation, Eq. dl It . was used in the last 
line. This proves that if Eq. d 1 3I > is valid for n, it is also 
valid for n + 1, which completes the induction step. Now, no 
particular use was made of the fact that $o stands for the black 
bead. We can as well start with bead n + 1 and work our way 
back to bead 0. In this way we obtain: 

$„ +1 = f/ n $l - t/ n _!$ . (15) 
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A fundamental, special case of this equation relates the am- 
plitudes of neighboring black beads connected by a spacer: 

$ s+ i - U s <f>! - E7._i$o. (16) 

Let us now turn to a dangling chain of k beads, whose ampli- 
tudes are $1, . . . , Now, from Eq. (|8) one has 

(2a;-l)* fc = (l-A)$ fc = $ fc _i. (17) 
Moreover, from Eq. (1151 we have 

*fe = £4-1*1 - tffc-2*0 (18) 

and 

$fc-l = C4- 2 $i - U k - 3 $ . (19) 
Inserting these into Eq. dl7l we obtain: 

(2x - l)(l/ fc _i$i - C4- 2 $o) - l/fc_2*i - C4-3$o (20) 
and hence: 

[(2a;-l)l7 fc _i-^ fc _2]$i = [(2z-l)£/ fe _ 2 -£/ fe _ 3 ]$o. (21) 
Making now use of Eq. ( II It it follows 

(C4 - I7*_i)$i = (C/fe-i - C/fc- 2 )*o. (22) 

Multiplying this equation with + Uk—i) and making use 
of Eqs. (IA6> and (IA7> of Appendix A leads to our second 
fundamental relation: 

U 2k *i = {U 2k -i + l)*o- (23) 

Now, let us consider a general black bead of amplitude $o 
that is connected to fx neighboring black beads (with ampli- 
tudes 3>a+i' • • • ) ^s+i) v ^ a cna i ns of s — 2k beads and that 
has attached to it f 2 dangling chains of k beads each, see Fig. 
[2 Then the black bead is /-coordinated, f = fx + f 2 - Equa- 
tion © corresponding to this situation reads: 

(h + h <tA)$o = J2 $ i Q) + E $ . (24) 

a=l b=l 

where, in an obvious way, $^ denotes the amplitude of the 
first bead of chain a. We now multiply Eq. i24\ by U s and 
insert the fundamental Eqs. (1161 and (1231 into it, obtaining 

h 

a=l 

+ h{U,-x + l)*a. (25) 

This can be rewritten as 

h 

[/ 1 -P/,,(A)]$ = E $ i+ ) i' ( 26 > 

a=l 



with 

P flS (\) = f + fU s - 1 -(f-a\)U s . (27) 

Here we recall that C/ s is a polynomial of order s in the vari- 
able x = 1 — A/2 (and analogously for {7 s _i), so that Pf tS (X) 
is a polynomial of order s + 1 in A. Equation i26\ is the funda- 
mental result for our model and it relates directly neighboring 
black beads; the white beads have been eliminated from it. 
We discuss it extensively in the next subsection. Evidently, in 
general one assumes that U s ^ 0, otherwise the multiplica- 
tions in Eqs. d23l and d25> lead to a loss of information. In 
the following we will take particular care of the cases where 
U. = 0. 



B. Determination of the eigenvalue spectrum 

Let us interpret this remarkable result, Eq. (I26> . This equa- 
tion relates the amplitude of every black bead to the ampli- 
tudes of its neighboring black beads only; all reference to the 
white beads on the spacers and on the dangling chains has 
been eliminated. We denote the eigenvalues of the reduced 
black bead system by A;, and have, rewriting Eqs. fl26i and 

h 

(/i-A 6 )<I»o = E*2i (28) 

a=l 

with 

Pf,.(A) = A b . (29) 

Remarkably thus (as in the case of the iterative procedure used 
for fractal nets 171 12110 each of the N — 1 nonzero eigenvalues 
Ab 7^ of the black subsystem leads via Eq. ( I29> to s + 1 = 
2k+ 1 eigenvalues of the full network (white and black beads). 
These eigenvalues correspond to eigenmodes of the system 
that involve the concomitant motion of white and black beads. 

In this way we have determined, given that our cluster is 
composed of N stars, a total of 

N P = (N- l)(2fc+ 1) (30) 

eigenvalues. We recall that the whole system must have 
Nfuii = N(fk + 1) eigenvalues, see Eq. (0}. Evidently, 
the remaining (in general, many) eigenvalues are connected to 
the zero eigenvalue, Af, = 0, of the black subsystem. 

The question arises to determine all the modes of the full 
system of white and black beads which lead to Ah = 0. Now 
in this case, as is known from general considerations M il Oil , 
the amplitude of all black beads is the same. This fact is also 
reflected by Eq.(l28>, which clearly admits as a solution 

^4 = 40. for a=l,...,h. (31) 

The simplest mode with Ab = and $o 7^ is the translation 
of all the beads in parallel; it corresponds to the eigenvalue 
Ai = of the full system. However, A& = is also related 



5 



to other modes in which the black beads move in parallel, and 
gives rise to additional eigenvalues of the full system. 

We now turn to consider the black subsystem. The simplest 
situation occurs when all black beads are at rest. In order to 
avoid a trivial solution one has imperatively to ask that (at 
least) one of the white beads neighboring a black one has as 
amplitude $1 / 0. Since the amplitude of its neighboring 
black bead vanishes, $o = 0, this implies, as a necessary 
condition for non-triviality, that U s = U%k = holds both 
for the spacers, Eq. d!6t (since also & s +i = 0), as well as 
for the dangling chains, Eq. ( 1231 . We are hence in both cases 
having U s = 0; recalling the form of U s given in Eq. iA7\ of 
Appendix A this means: 



U s = U 2k = (U k + U k -i){U k - U k -x) = 0. 



From the point of view of the spacers, this implies that 
the eigenvalues corresponding to Eq. fl38i have eigenmodes 
which are antisymmetric with respect to the mid-chain beads 
(k, k + 1), whereas the eigenvalues corresponding to Eq. d37l > 
have symmetric eigenmodes. 

The modes which we have now determined and which cor- 
respond to U s — are highly degenerate. We display the 
evaluation of their (corresponding) multiplicities in Appendix 
B. The basic idea developed in Appendix B is that, due to Eq. 
(|6j, for $o = 0, a non-vanishing value of $^ a ' cannot appear 
alone, but must be compensated by the motion of (at least) 
another bead. Hence, if one starts to look for linearly inde- 
pendent eigenmodes it is reasonable to normalize one of the 

(32) amplitudes, and set $^ = 1. Moreover, one has to satisfy 



Moreover, consideration of Eq. J22I shows that for dangling 
chains one has to sharpen the condition and to ask even that 



U k - U k -! 



0. 



(33) 



From Eq. ( 1321 we thus obtain two kinds of normal modes. 
One set of them obeys Eq. J33i and may involve both dan- 
gling chains and also spacers, whereas the second set of nor- 
mal modes obeys 



U k + U k -! 



0. 



(34) 



and can be fulfilled only by modes localized on the spacers 
exclusively. 

Now the roots of U s = U 2k = are, using Eq. d!2i . given 
by sin[(2fc + l)ip] = under the condition that sintp ^ 0. 
Hence they have the form 



2k + 1 



with r = 1, . . . , 2k. 



(35) 



From this relation the roots x r follow from x r = cos<p r . 
Given the form of Eq. J35i it is clear that none of them equals 
unity, x r ^ 1, and thus none of them leads to A r = 0, see Eq. 
©• 

Moreover, not all of these x r (as stressed above) may ap- 
pear when white beads of dangling chains move. For these 
Eq. ( 1331 . namely 



U k (costp) - C4_i(cosip) = 



(36) 



must hold. From Eq. d!2i the roots of Eq. d36l are the roots 
of sin[(fc + l)(p] — sm(k(p) = for which sinip =^ 0. Since 
sin(a + 0) — sin(a — (3) = 2 sin (3 cos a, the corresponding 
roots are: 



Pm = 



2m -1 
2k +1 



with m — 1, . . . ,k, (37) 



whereas the modes involving spacers can also have, see Eq. 

C3, 



2m 
2k + l' 



with m = 1, . 



,k. 



(38) 



A small remark concerns now the symmetry of the modes con- 
sidered. From Eq. dl 3i we have, since $o = that 



Eq. 0. Again, due to simplicity, one assumes that all the 
other white neighbors of the given black bead are at rest, ex- 
cept for another one, or, at most two. This leads to having, 



in an obvious notation = —1 or &\ 0> = §\ c ' = -1/2, 
respectively. In this way one is able to count all the linearly 
independent pathways through the full system, in which the 
signs of the modes are alternating, and which, depending on 
the eigenvalue considered, may of may not involve dangling 
chains. 

This procedure is performed in Appendix B, where we 
show that the multiplicities of the eigenvalues to Eq. d37l > and 
to Eq. fl38i are (quite generally) given by: 



s(b) 



do) 



and by 



A+ = l-(N-l) 



(f-l)N-l, 



(40) 



(41) 



where the indices indicate that they belong to Eqs. d34l > and 
( I33l l. respectively. As before, in Eqs. d40l > and ( 14 11 N denotes 
the number of stars and / the number of additional links con- 
necting them. Then A + is related to the number of loops in 
the system. For A + = the structure is topologically a tree, 
no loops exist, and eigenvalues to Eq. d38l do not appear in the 
spectrum. For A + = 0, furthermore, A _ = (/ — 2)N +1, and 
A_ + 1 is the number of dangling chains. From Eqs. d40l > and 
( 14 H the total number of eigenvalues arising from this class of 
modes is 



iV A = fc(A+ + A_) = k[(f - 2)N + 1]. 



(42) 



As a side-remark, we note that in a very special situation 
(seldom encountered in randomly linked systems), for a par- 
ticular value of A;, in Eq. ( I29l > k modes of the black and white 
system get to be such as to have U s = 0. The corresponding 
modes belong then to both classes, Eq. d29i and Eq. ( I4U . 
(Evidently, they may only be counted once). The situation is 
discussed in Appendix C. 

Finally we have to consider the situation in which all black 
beads are moving in parallel and the motion of the white beads 
compensates it. Then usually U s ^ and one is again led to 
reconsider P/, S (A), but now, since = one has 



*n/C/n-l = ®n+l/U n 



(39) 



Pf,sW=0- 



(43) 



6 



Moreover, using the relation s = 2k and Eqs. (IA6t and JA7> 
of Appendix A with n = k, we remark that Pf a (A), Eq. 1271 . 
can be reformulated as 

Pf,.W=Qf.'W(U k + U k - 1 ), (44) 

where we set 

Qf.sW = f(U k -i - tffc-a) - (/ - ^)([/ fe - C/^O. (45) 

From Eq. (1441 and recalling the remarks after Eq. ( I27l i. 
P/ iS (A) = requires 

Q/>(A) = 0, (46) 

whose modes involve spacers and dangling chains alike. As a 
side-remark, it is sometimes convenient to formulate (3/, s (A) 
in an alternate way. From the recursion relation, Eq. il Q , one 
has 

(Uk^-Uk^-iUk-Uk-x) = 2U k . 1 -2xU k+1 = \U k - lt 

(47) 

with which Eq. ( 145 \ takes the form: 

Qf,,{X)=X[{f-a)U h - 1 + aU k ]. (48) 

Turning now to the general features of Q f tS (A) we remark that 
its degree, as given by Eq. J45b . is fc + 1, and hence we expect 
in general 

N Q = fc + 1 (49) 

different eigenvalues from Eq. J46t . By inspection of Eq. 
J48t it is evident that Ai = is indeed a solution of Eq. < !46b . 
Hence the translation of the full system, which is related to 
Ai = 0, is included in the set of modes corresponding to Eq. 
J46i . This solution, as well as the other fc additional ones of 
Eq. J46V is non-degenerate. 



linearly independent eigenmodes, where in the last expression 
we recalled Eq. ©. By this we have indeed found a complete 
set of eigenmodes of the full system, leading to the correct, 
total number of eigenvalues. 

To end this section we stress that the determination of the 
multiplicities Eqs. d41i and J40i for the modes given by Eqs. 
03t and 04l >. i.e. to the relation U s = 0, see Eq. d32i . was 
very general and depended only on the topology of our net- 
work. Eqs. J41i and d40l were discussed in the framework 
of normal modes in which all black beads are at rest; the de- 
generacy of their eigenvalues is inherent from the network's 
structure. 

Interestingly, in special cases (special types of lattices or 
parameter values), one may also be led to the relation U s = 0, 
albeit for the other two types of normal modes. One obtains 
then additional degeneracies, which are, however, acciden- 
tal. Thus for lattices without dangling chains and where the 
black beads build two alternating sublattices (see Appendix 
C) eigenvalues to S (A) = Ab = 2/ turn out to be of the 
type U s — 0. Furthermore, the choice a — f/2 for the pa- 
rameter a transforms Eq. ( I46> . Q/,s(A) = 0, using Eq. J48l > 
to the condition Uk-i + Vh = 0, which automatically implies 
U s = U 2 k = 0, see Eq. ( IA71 . 



IV. APPLICATIONS 

In this Section we apply our method to several types of lat- 
tices made out of star polymers. In order to exemplify our 
ideas we start with the simplest case, namely with a single 
star polymer. 



A. Star polymer 



C. Classification of eigenvalues and special cases 

Summarizing, we have found a set of equations that fully 
determine the eigenvalue spectrum of our system. These are, 
namely, Eqs. (j29j, J33l >. J34b and J46I . Evidently, they belong 
to different types of modes, which are clearly linearly inde- 
pendent of each other. In the first class, given by Eq. ( I29l >. 
the coordinates of some of the black beads are non-vanishing 
and the black beads do not move in parallel. The second class, 
given by Eq. ( 146ft . is that in which all the black beads move in 
parallel. The next two classes, exemplified by Eqs. J33l > and 
04> . are such that all the black beads are at rest. 

Let us check the total number of eigenvalues that we get 
from these modes. We find from Eqs. ( I29l l. ( 142b and ( I46i a 
total of 

N P + N Q + N A = 

= (2k + 1)(N - 1) + (fc + 1) + k[(f - 2)N + 1] 
= (/fc + 1)N = Nf u u (50) 



As a first example of our formalism we consider a star con- 
sisting of a black core out of which emerge / branches of fc 
white beads each. Evidently, in this case the black system 
consists of a single bead, with the unique eigenvalue A5 = 0; 
then Pf tS (X), Eq. ( I29K does not appear. Furthermore, from 
Eq. ( 1401 1 one has A + = 0, by inserting N = 1 and I = 
into it. Apart from Ai = there are two classes of non- 
trivial solutions: one class stems from Q/ S (A) and in it the 
black bead moves (the modes are called symmetric in the lit- 
erature l22Tl . see Refs. fT^ l23ll for the situation for den- 
drimers); in the other class, stemming from the A_ case, the 
black bead is at rest (this leads to antisymmetric modes 1 22]). 
From our formalism (and as it is well-known) the symmet- 
ric modes are non-degenerate; there are fc of them and they 
fulfill Eq. i46i . The antisymmetric modes are, on the other 
hand, A_ = (/ — 1) - fold degenerate; they lead each to the 
fc distinct eigenvalues given by Eq. d37t . For the star the rea- 
son for the degeneracy of the antisymmetric modes is evident; 
it can be seen by pairwise combining in linearly independent 
manner, the / branches of the star. 
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B. The Cube polymer 

Another very interesting example for the use of our formal- 
ism in a simple context is provided by a connected structure 
topologically identical to a cube. The basic network (black 
system) consists here of eight black beads which are the cores 
of eight stars with / = 3 arms each and with k white beads 
on each arm. The arms are then connected pairwise, through 
I = 12 links, see Fig. 0. This leads to 12 spacers of s = 2k 
white beads along the "edges" of the "cube". 




FIG. 2: The Cube polymer. As a general remark on its structure in 
solution, see caption of FigQ 

The cube contains loops of spacers but no dangling bonds. 
Moreover, the structure is alternating, in the sense of Ap- 
pendix C. Hence, a very interesting aspect of the cube is that 
it provides a very nice example of the special case treated in 
Appendix C: A rare, additional degeneracy related to Pf tS (A) 
and to U s = occurs. 

For the cube the eigenvalues of the black system are 6,4,2 
and 0. Of these, A& = and A& = 6 are non-degenerate, 
whereas A& = 4 and Af, = 2 are each three-fold degenerate. 
Thus, our formalism leads, based on Eq. 1291 . to 2k + 1 solu- 
tions to Pf tS (A) = 6 and to 3(2fc+ 1) solutions to P/, s (A) = 4 
and Pf )S (A) = 2, each. We also find from Eq. J46> k + 1 
solutions corresponding to the equation Q/, s ( A) = 0. Fur- 
thermore, from Sec. 3.2 we compute, based on Eqs. d40> 
and J41> . the multiplicities of the k modes to Eqs. ( I34> and 
(|33} to be A+ = 5 and A_ = 4. All in all, we thus obtain 
7(2k + 1) + (k + 1) + 9fc = 8(3fc + 1) modes. 

We remark now that for the cube one can construct a loop 
on every of the six faces, but that from these six loops, there 
are only five which are linearly independent. As far as each of 
these loops has an even number of chains, it can accommodate 
modes corresponding to Eq. (I37> that are symmetric at mid- 
spacer. Thus, the total multiplicity of these modes is A_ + 1 = 
5. As shown in Appendix C, the additional multiplicity of 
these k modes is already accounted for by k of the solutions 
to Pf_ s (X) = 2/ = 6. For these k modes, the factor U2k 
vanishes, see Appendix C, meaning that the black beads are at 
rest. We recall that this is a general feature of fully connected 
networks, in which all loops have an even number of chains, 



see Appendix B. 



C. Periodic networks 

As a third example we now apply our formalism to a topo- 
logically regular d-dimensional network, as it is treated in de- 
tail for d = 3 in Ref. Q. 

We start with a periodic, lattice-like arrangement of N — 
x x • • • x /yW black beads and impose periodic 
boundary conditions. By this all the black nodes get to be 2d- 
coordinated, / = 2d. The topological regularity of the system 
of black beads allows readily to determine 1 24] its eigenvalues 
At,. One finds: 

A 6 = 2d(l-cose 6 ), (51) 
where cos 0t,, the so-called generating function is 

1 d 

cos 6b = — costfj, (52) 
d 1=1 

with <&i = Si(2n)/N^ and Si = 0,..., N (i) - 1. One may 
note that Eq. {51} provides all the x x • • • x 
eigenvalues of the black system and that the zero eigenvalue, 
Ai = 0, corresponds to the choice Si = S2 = ■ ■ ■ = Sd = 0. 

We now replace all bonds of the black network by spac- 
ers with 2k white beads each. At this stage one has to solve 
Eq. d29t to obtain the first class of modes of the whole (black 
and white) network. Now, as has been noted in Ref. 12411 . 
a considerable simplification of the problem occurs when the 
coefficient a happens to equal //2, i.e. for a = f/2 = d. 
One may see the idea clearly by starting from Eq. (1271 for this 
parameter choice and recalling Eq. (|9j: 

P f , s (X) = /[i -(1-^)17, + [/,_!] 

= f(l-xU„ + ?7 s _i). (53) 

The condition Eq. d29l is then, with Eq. ( I5U : 

cos6& = xU s — U s -x- (54) 

Now, using Eq. d!2i . it follows 

cos 8b = {cos(<^) sin[(s + l)ip] — sin (sip)}/siriip = 

= cos[(s + l)^, (55) 

since 

sin (sip) — sin[(s + X)ip — </?] = 

sin[(s + l)ip] cos ip — cos[(s + l)ip] sin ip. (56) 

Equation d55i is readily solved, by which we recover the so- 
lutions of Ref. |j2l: 

^m = ^77±-TT' with m = 0, 1, 2, . . . , fc, (57) 

s + 1 s + 1 
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FIG. 3: A second generation dendrimer with spacers. Here, g = 
2, f = 3 and each spacer chain contains two beads (s — 2,k — 
1). The dash-dotted line indicates one of the stars from which the 
structure is built. As a general remark on the dendrimer's structure 
in solution, see caption of Fig.Q 



where one has to keep only one solution for m = 0. Excluding 
the case X b = 0, Eqs. (|52} and {55} lead to N P = (N— l)(s+ 
1) eigenvalues. 

We now turn to the degenerate eigenvalues corresponding 
to Eqs. (37} and (|38}. Using Eqs. (glj and (gl} for A+ 
and A_ and recalling that we have here a highly connected 
periodic network, without dangling bonds, whose number of 
links is I = dN, we find as in Ref. Q A_ = (d - 1)N and 
A+ = A_ + 1. 

It remains to study Eq. J46i . Now for the choice a = f/2 
one again has with Eq. J48l > the simplification 



Q f , s (\) = d\(U k + U k - 1 ), 



(58) 



which leads to an additional, accidental increase of the num- 
ber of eigenvalues corresponding to Eq. (1381 . 

Finally we note that our approach is not limited to peri- 
odic boundary conditions, but can handle as well finite regular 
structures with dangling half chains. 



V. DENDRIMERS WITH SPACERS 

As a specific example of branched polymers, we treat here 
star burst dendrimers, which are built from /-functional stars, 
linked in such a way as to describe a regular Cayley tree, see 
Figure[3] One may remark that our renormalization procedure 
reduces the dendrimer with spacers and dangling bonds of Fig. 
|3]to a simple Cayley tree. Now, let us recall the properties of 
the eig envalue spectrum of such a simple dendrime r (T^L fl8l 
[Hill IH El- As has also been noted in Ref. (iK this 
problem can be partly mapped to a rescaled one-dimensional 
Gaussian chain. Thus, it allows for a formulation in terms of 
Chebyshev polynomials 1 15]. 



We start with the system of black beads, i.e. with the den- 
drimer at the generation g. The total number of black beads is 
then: 



N = 1 + f 



if - l) g - 1 
f-2 ■ 



(59) 



As discussed in Refs. 1 16:. ll8lll9U23ll the eigenvalue spectrum 
corresponds to two different types of modes: those in which 
the core is mobile and those in which the core does not move. 

Equations for the eigenvalues of the classical dendrimer 
of generation g and general functionality / have been derived 
in Refs. 1191 12311 . In terms of the Chebyshev polynomials of 
the second kind these eigenvalues can be obtained from the 
roots of the following equations, see Appendix D: 



U g (x) + 



(60) 



and 



U m (x) - v7- lU m -i(x) = 0, with m = 1, . . . ,g. (61) 



The eigenvalues themselves are then given by 



A 6 = / - 2xy/f^l. 



(62) 



We remark that for g + 1 < v\7 — 1 and for m < (to — 
1)V / — 1 these solutions include the so-called exponential 
modes with |a;| > 1, for which Eq. (1 1 2i may be written in 
terms of hyperbolic functions. The g solutions to Eq. (I60> 
are non-degenerate, while the degeneracy of the eigenvalues 
to Eq. ( 16 U is given by 



A, 



(/ ~ 1) if to = g 

/(/-2)(/- otherwise _ 



(63) 



This defines the full eigenvalue spectrum of the dendrimer 
consisting of the black beads only. 

Now, we calculate the eigenvalues of the full system (the 
system including spacers and dangling chains) by solving Eqs. 
j^9l , d3H and gS). Each eigenvalue X b ^ of the black 
system gives rise to a set of s+ 1 eigenvalues of the full system 
which inherit its degeneracy. The degeneracy of the solutions 
to Eq. (1331 is A_, see Eq. (14 1 1 , while the solutions to Eq. 
d46l are non-degenerate. 

This two-step approach (of first finding the eigenvalue spec- 
trum of the reduced black system, and only then calculating 
the eigenvalue spectrum of the full system, by solving Eqs. 
( I29> , d33l and (I46», considerably reduces the complexity of 
the problem. We should note that the maximal degrees of 
the polynomial equations involved here are g and s + 1, As 
shown below, this allows us to numerically treat dendrimers 
with spacers having a total of more than 14000 beads. This 
demonstrates one of the main advantages of our present work, 
when comparing it to a one-step approach 1 15]; such an ap- 
proach leads here to a polynomial equation of much higher 
degree, namely of degree (s + l)(g + 1). We also note that 
in Ref. [15] the peripheral beads of the dendritic system were 
assumed to be fixed in space. In retrospect, this assumption, 
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FIG. 4: Storage modulus G'(uj) for a dendrimer of generation g — 6, 
functionality / = 3, with spacer chains of length s = 4 in between 
the branching points and dangling bonds of length k — 2. Here, a 
varies from 0.01 to 100. 
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FIG. 5: Loss modulus G"(ui) for a dendrimer of generation g = 6, 
functionality / = 3, with spacer chains of length s = 4 in between 
the branching points and dangling bonds of length k = 2. Here, <j 
varies from 0.01 to 100. 



In Figs. [4]and[5]we display G'{u>) and G"(u>) as functions 
of the reduced frequency Q for dendrimers with chain length 
parameter s = 2k = 4 while varying cr from 0.01 to 100. 
The mobility of the white beads (i.e. Co) is kept constant. The 
curve corresponding to cr = 1 describes the case when all the 
beads in the system have the same mobility. As it was shown 
frU \T% . differences in the mobilities of the different kinds 
of beads lead to the appearance of plateau-type behaviors in 
G'(u>) and to additional peaks in G"(u>). For a > 1 the spac- 
ers and dangling beads are the "lighter" beads, being more 
mobile than the branching units (black beads). Their num- 
ber being six times larger, their relaxation dominates the spec- 
trum. In particular, for a = 100, the loss modulus G"(u>) dis- 
plays a major peak at intermediate frequencies, corresponding 
to the relaxation of the white beads and a second minor peak 
at lower frequencies, corresponding to that of the black beads. 
For a < 1 the situation is reversed: The black beads are more 
mobile and the minor peak corresponding to their relaxation 
is shifted to higher frequencies, while the major peak induced 
by the white beads keeps its position. Corresponding effects 
are seen also in the behavior of the storage moduli G'(uS), 
which for a ^ 1 display a plateau resulting from a gap in 
the relaxation. The plateau is formed between the domains 
of high and of low frequency relaxation. For a = 100 the 
white beads dominate the high frequency relaxation, leading 
to a lower level of the plateau than for a = 0.01, where the 
high frequency relaxation of the black beads is only weakly 
evident. 
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besides of being rather arbitrary, appears to be counterproduc- 
tive. 

Let us now demonstrate numerically the effects of varying 
the parameters s = 2k and cr, which determine the length of 
the spacers and of the dangling chains, and the ratio of the 
friction coefficients, respectively. We choose here a g = 6 
generation dendrimer with functionality / = 3, and show its 
storage G'(u>) and loss G"(oj) moduli, calculated according 
to Eqs. (|2ji and 0. 



FIG. 6: The storage modulus for a dendrimer of generation g — 
6, when the number of spacers between the branching points grows 
from s — 2 to s = 50. In all the situations a — 100. 

To demonstrate the influence of increasing the length of the 
spacers and of the dangling chains, we plot in Figs. [6]and0the 
corresponding results for the moduli G' (u>) and G" (u>). Again 
we consider a g — 6 generation dendrimer with functionality 
/ = 3. Fixing the parameter a = 100 (higher mobility of 
the white beads) we vary the length of the spacer and of the 
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dangling chains, by increasing s — 2k from s = to s = 50. 
For s = 50, the full system consists of Nf u u = 14440 beads, 
while the reduced (black) system is a dendrimer containing 
N = 190 beads. The classical dendrimer (s = 0) shows 
a typical, non-scaling behavior of the relaxation part of the 
spectrum, both for G'{uS) and for G"{lo). Increasing s the 
relaxation moduli display a mixture of dendrimer and chain 
like behavior. For short chains we observe in G"{lo) a sec- 
ond minor peak at lower frequencies, which, for longer chain 
lengths becomes less and less pronounced, remaining eventu- 
ally only visible as a shoulder. For long chains both G'(ui) and 
G"{uj) are dominated by the chain like relaxation behavior of 
the white beads. 




FIG. 7: The loss modulus for a dendrimer of generation g — 6, when 
the number of spacers between the branching points grows from s = 
2 to s = 50. In all the situations a — 100. 



VI. CONCLUSIONS 

The dynamics of many branched polymer systems can be 
understood in terms of the dynamics of end-linked star poly- 
mers. Here, we presented a general method to analytically 
treat this problem, while also including the important effects 
brought along by changes in the length and in the mobility of 
the arms of the stars. We developed a mathematically exact 
procedure which allows to determine the full eigenvalue spec- 
trum of such networks, by exactly mapping the eigenmode 
spectrum of a network formed by end-linking star polymers 
to that of a reduced, core system. Namely, for any such net- 
work (under the weak restriction that adjacent stars are linked 
by one bond only) one obtains the core system by replacing 
each star polymer by its core. Given the eigenvalue spectrum 
of this core system we obtain the full spectrum of eigenval- 
ues by solving a set of polynomial equations. For a network 
of star polymers with / arms of k beads each, every non-zero 
eigenvalue of the core system gives rise to 2k + 1 non-trivial 



eigenvalues of the full network. These constitute the first class 
of eigenvalues; the second class is related to the zero eigen- 
value of the core system and is organized in three groups: The 
first two groups correspond to the symmetric and antisymmet- 
ric modes of the arms of a single star polymer, while the third 
group corresponds to the modes of spacer chains with fixed 
ends. For all these modes we determined their eigenvalue 
multiplicities; these depends only on the topology of the core 
system, i.e. on the total number of stars and on the number of 
independent loops in the core system, as worked out in Ap- 
pendix B. 

We demonstrated the full power of our approach by recov- 
ering readily the full spectra of model structures such as star 
polymers and regular <i-dimensional networks; we reproduced 
former results for d = 3, and generalized these to arbitrary d. 
Furthermore, we exemplified hidden symmetries of the prob- 
lem using the "cube" polymer as fundamental example. Ad- 
ditionally, we performed extensive calculations on dendritic 
structures containing spacers and dangling chains. Our two- 
step approach allowed us to readily treat systems containing 
more than 14000 beads; in fact, numerically we are limited 
only by our possibility of determining the spectra of the core 
system; this means for arbitrary systems, of some 5000 black 
beads, and (again) practically unlimited for hierarchical (frac- 
tal, dendritic, regular) core networks. We note that a previous 
approach to dendrimers with spacers lH5ll . restricted to sys- 
tems with fixed peripheral junctions, did not succeed in de- 
coupling the problem and obtained only approximate spectra 
(sometimes, even negative, unphysical eigenvalues). 

Considering the relaxation moduli of copolymers, we find 
strong effects when the mobilities of the beads involved differ. 
Also, increasing the length of the arms of the stars leads to 
smooth transitions between the relaxation behavior due to the 
core and due to the linear chains. 

Our method is very general, in that the topology of the core 
and the details of the stars' geometries are (in a large measure) 
arbitrary. It provides a scheme which considerably reduces 
the calculational efforts, by splitting the task into two steps. 
Furthermore, if any symmetry is present in the core system, 
additional reduction steps are possible. As an outlook, we 
note that the method presented here allows to generalize re- 
cent analytic results for the dynamics of randomly branched 
polymers ll27l EH E9I1 to systems built from them by the in- 
sertion of spacers and of dangling bonds, and in which the 
mobility of the beads may vary. 
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APPENDIX A 

In this Appendix we derive several relations connecting the 
Chebyshev polynomials of the second kind, U n (x), Eqs. (lilt 
and d!2i of the main text. We start from the relation: 

sin (2mp) sin cp = 2 sin (mp) cos (mp) sin ip = 

= sin (rup) [ sin (mp + ip) — sin (n</9 — ip) ] 

= sin (n< / 9){sin \(n + l)ip] — sin \(n — l)p]} (Al) 

Using Eq. (I12L it follows by division with sin 2 i/p that 



Uzn-i — U n -i(U n — U n - 2 )- 



(A2) 



Moreover from: 



sin [(n + l)ip] sin [(n — l)ip] = 

= [cos ip sin (mp) + cos(mp) sin ip] x 

x [cos (/? sm(nip) — cos(n<p) sin </?] = 

= [cos 2 ip sin 2 (nip) — cos 2 (ny>) sin ip] 

— cos 2 99 sin 2 (nip) — [1 — sin 2 (nc/?)](l — cos 2 (/?) 

= cos 2 tp sin 2 (ny) — sin 2 ip + sin 2 (mp) — cos 2 ip sin 2 (mp) 

= sin 2 (nip) — sin 2 <p, (A3) 

it follows by division with sin 2 ip that 

U n Vn-2 = Ul-Y - 1, (A4) 

and hence 

C/ 2 _! - [/„[4-2 = 1. (A5) 
Thus, from Eqs. JA2t and ( IA5> 



(tf n + [/ n _!)(C/ n _i - tf„_ 2 ) = U n U n -! + U 2 _, - 

-U n U n - 2 - U n - X U n -2 = U 2n -i + 1. (A6) 

Furthermore, one also has 

(U n + U n -i)(U n -U n - 1 ) - U 2 - U 2 _ ! = £/ 2 n , (A7) 

as we proceed to show based on the relation 

sin 2 991 — sin 2 ip 2 = sin (t/?i + p 2 ) sin (tpi — ip 2 ), 

Eq. (4.3.20) of Ref.EjJl. Thus 

sin 2 [(n + l)ip] — sin 2 (nip) — sin[(2n + l)ip] sin tp. (A8) 

Dividing by sin 2 <p and applying Eq. d!2i proves now Eq. 
dA7}. 

APPENDIX B 

In this Appendix we derive the multiplicities, Eqs. (I40i 
and (14- 1 i of the modes to Eqs. J33i and (I34> . by an explicit 
construction. Let $ and ^ be modes on a spacer with 2k white 
beads with fixed end black beads, to given eigenvalues, Eqs. 
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FIG. 8: Directed paths Vi,...,V5 between the dangling chain 
(M + 1) = 6 and m (m = 1, . . . , 5) and the corresponding modes 

$( m J = ±$ for a structure obtained by end-linking N = 4 poly- 

I a 

mer stars. Indicated by +/— signs are the amplitudes of the modes 
<2>( m ). 



d3 3I > and ( I34K where $ is symmetric and ^ is antisymmetric 
at mid-spacer: 



fc+ij 



and 



fe+1- 



(Bl) 



Consider now one full structure (with white and black beads), 
formed from end-linking stars by I links, I > li = N — 1. 
If li — I then the structure is topologically a tree and the total 
number M + 1 of dangling chains is M + 1 = (/ - 2)N + 2. 
For I > 1 1 the network contains loops and we set L = I — h, 
so that I = L + JV — 1. By removing now from the struc- 
ture L judiciously chosen links out of the I links, the ensu- 
ing graph turns out to be a connected tree. We number then 
the new dangling chains (created by removing the L links) by 
(1, 2), (3, 4), . . . , (2L - 1, 2L) and the other dangling chains 
of the tree (these exist for 2L < M- 1) by 2L+1, . . . , M+l. 

We now assume 2L < M (then M + 1 is a dangling chain 
of the original structure) and specify M paths Vm, where V m 
is the directed path between the dangling chains M + 1 and 
to (with to = 1, . . . , M), Since our structure is a tree, each 
directed path is unique. For each path V m we define a mode 
<j>( m ) by its values when restricted to any (half-) spacer a £ 



$0) 



(B2) 



where the sign oscillates between adjacent spacers along the 
path (see Figure |8j. Here we use an implicit numbering (in 
the direction of the path V m ) of the beads along each spacer, 
with obvious restrictions for the dangling chains. We fix the 
+ sign for the dangling chain to. Outside the path Vm, the 
mode $( m ) vanishes. Each mode $( m ) is formed by symmet- 
ric pieces along the spacer and fulfills Eq. J24i at the black 
beads. The modes $( m ) are obviously linearly independent, 
since each path V m contains at least one segment that is not 
contained in any other path. In the case A + =0we have thus 
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constructed M — (f ~ 2)N + 1 different degenerate modes 
corresponding to a given eigenvalue of Eq. A37I . by which we 
have proven the correctness of Eqs. J40i and d41t in the case 
of a tree structure. 

In the following we start to reconnect the dangling chains 
(1, 2), (3,4), etc. and follow the evolution of the system of 
eigenmodes under this operation. For this we define for each 
path V m on the tree a mode vf^" 1 ' by its restrictions 



(m) 



with a € V n 



(B3) 



while we let it vanish outside V m . Note, however, that, due 
to its eigenvalue, <3>( m ) cannot involve the motion of dangling 
chains, a fact which we ignore for the moment. Furthermore, 
we construct the modes: 



<jj(m,m+l) _ <j)(m) _|_ <jj(m+l) 



and 



(m,m+l) \T;( m ) 



(m+1) 



(B4) 



(B5) 



with to = 1, 3, . . . , 2L — 1 The union of paths V m U V m +i, 
consists of a loop C m and a "tail" T m between this loop and 
the dangling chain M+1 given by the intersection V m C\V m +\ 
(see Figure [9jl. We will call C m an even or odd loop, if it 
contains an even or odd number of spacers. 




m=l 



(M+l)=6 




FIG. 9: The loop £3 and tail T3 obtained when introducing one ad- 
ditional link between the dangling chains m = 3 and m — 4. We 
also indicate the amplitudes of the mode $' 3,4 ^ by +/— signs. 

By construction, the modes ^(jn,m+i) and fr(m,m+i) f u i_ 
fill the boundary condition JB U at the center of the spacer 
and Eq. J24> at the branching point of paths V m and V m +i- 
Mode ^r( m ' m + 1 ) is non-zero only on the loop C m , whereas 
mode $v m > TO + 1 ) is also non-zero on the tail T m if C m is odd. 
Note that each linkage (reconnection) reduces the number of 
<f>-type eigenmodes by unity and increases the number of it- 
type eigenmodes by unity. By insertion of L links one obtains 
L modes belonging to Eq. JB4L L modes belonging to Eq. 
iB5\ . but has now only M — 2L modes $( m ) on the paths ? ra 
(m > 2L) connecting dangling chains. We obtain thus: 



and 



A 1 



(M — 2L) + L = M — L. 



(B6) 



(B7) 



Given that M = (/ - 2)N + 1 and I = L + N - 1, these 
expressions are equivalent to Eqs. j40t and J41l i. 

We turn now to the case which we excluded up-to-now, 
namely to 2L = M + 1. Then the original structure has no 
dangling chains. In this case we perform first the steps above, 
and reconnect (1,2), (3, 4) ... up to (M - 2, M - 1). By this 
only the path Vm between the dangling chains M and M + 1 
remains to be closed. Connecting these two dangling chains, 
the path Vm becomes a loop, which we denote by Cm- 
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FIG. 10: Two examples of situations encountered while constructed 
the modes w hen closing the last link (M, M + 1). 



First, we assume that some of the loops Cj for modes corre- 
sponding to Eq. JB5I have an odd number of spacers (Figure 
1 1 0i . Let the paths be ordered such that Cm is odd. The above 
considerations remain valid for all modes ty( m ' m + 1 ) as well 
as for all modes $( m > m + 1 ) with an even C m . However, upon 
closing the last link transforming the path Vm to a loop Cm, 
all modes that include a nonvanishing tail part Tj, i.e. 

those for which Cj is odd, have to be redefined. The last link 
consists of the two dangling half chains M and M + 1. The 
amplitudes of a given mode on the chains along the 

tail Tj are ±2: hence on the dangling chain M + 1 the mode 
takes the value 2ej, with either €j — +1 or ej = — 1. On the 
dangling chain M the amplitude of vanishes; closing 

the loop C m we are faced with an amplitude mismatch of 2cj 
between the dangling chains M and M + 1, fact which we 
have to mend. Now, since here Cm is an odd loop, the mode 
cj)(-M) (i e fi nec i on -p M nas amplitude +1 on the dangling chain 
M and amplitude —1 on the dangling chain M + 1. Thus, 
closing the loop Cm, the mode $( M ) will exhibit an ampli- 
tude mismatch of 6 = 2 between the dangling chains M and 
M+1, while being well defined everywhere else. This fact 
allows us to construct for each odd loop Cj a well defined $- 
type mode (which lives on the union of Cj, Tj and Cm) by 
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consider a general site i and its /j neighbors 



J 1 1 



(B8) 



where 2ej is the amplitude of mode $Cm+ 1 ) on thf dangling 
chain M + 1. For all j < M with an odd loop Cj these modes 
fulfill Eq. ( 1241 and all other boundary conditions. On the other 
hand, <J/( M ) becomes a valid mode on the loop Cm- Thus, as 
before, the closure of a loop leads here to the disappearance of 
one $-type mode $( M ) (counted by A_) and to the creation 
of a v^-type mode (counted by A+). Hence Eqs. dB6l and 
dB7l . and thus Eqs. go) and gD hold - 

Finally, we concentrate on the special situation when all 
loops C m are even. Then, when closing the last link, all modes 
dB5i remain valid, since no nonvanishing tail parts have to be 
considered. Formally, the closure of the loop allows now to 
define two new modes, see Eqs. iB4\ and (IB5> . one of $- and 
one of 'I'-type, which replace the previous $( M )-mode. How- 
ever, as shown in Appendix C, the additional <f>-type mode 
which would increase A_ by unity and would lead to k dis- 
tinct modes of the whole, black and white system, is only ap- 
parently new. All these seemingly new k modes are in fact 
already included in the solutions of Eq. ( 1291 . as we prove in 
Appendix C. 



APPENDIX C 

In this Appendix we study whether for certain At, ^ the 
relation Pf tS (X) = Xb may have solutions related to U s = 0. 
Here we show that this can occur only if Xb = 2f. Now, from 
Eq. d!2i . U s (cosip r ) = means that 



fr 



with r = 1,2, . . .,2k, 



2k + 1 
from which follows that 

sin [(2k + l)tp r — if 



(CI) 



[/ s _i(cOS tf r ) 



sin ip r 

sin (rir — ip r ) _ ( +1 for r odd 



sin ip r 



- 1 for r even 



(C2) 



For U s = 0, U s -i takes one of these two possible values. 
Making now use of Eqs. 1271 and d29t we see that: 



h = f+fU a -i 



2/ for r = 2m — 1, with m = 1, . . . , k 
for r — 2m, with m = 1, . . . , k. 

(C3) 

Hence, the only possibility is to take Xb — 2f, given that we 
are looking for solutions to P/ jS (A) = Xb with At 7^ 0. 

To show now that indeed A& = 2/ is a possible eigenvalue 
in the system consisting of the black beads only, we analyze 
the existence and the properties of an eigenmode to it. We 
start with Eq. J26l > for the system of black beads, where we 



(/i-2/)^ = 



(C4) 



(C5) 



a=l 



We now pick the site j such that on it |$j | is maximal, |$j| = 
maxi\<&i\. Furthermore we choose the amplitudes such that 
&j is positive. Then from Eq. dC5t 



2/$ 



fi 

a=l 



(C6) 



This means that, on one side 2/ < 2fj, whereas, by construc- 
tion, fj < f. It follows that we must have fj = f. 

Moreover, inserting fj = f into Eq. iCbt it also follows 
that we must have |$j — <f> jQ | = 2|$j| for all a, which sim- 
ply implies §j a — — $j for all the neighbors of j. Obviously, 
this argument can be repeated for all the / neighbors of j, etc. 
Hence, the only mode compatible with A^ = 2/ is an alternat- 
ing mode of two sublattices moving against each other. Such 
a mode is only possible when: first, all black beads have the 
functionality / on the black lattice, and secondly, all loops in 
the structure are even. One should note that the first condition 
means that no dangling chains are allowed. 

We now turn to the motion of the full lattice (white and 
black beads); we recover it by replacing each bond of the 
black lattice with a spacer of 2k white beads. From Eq. 
dC3l > the condition A^ = 2/ implies that the r must be odd, 
r = 2m — 1, and that there are exactly k solutions satisfying 
this condition. Furthermore, based on Eq. fl39i it follows that 
all these solutions correspond to eigenmodes which are even. 
Symmetric (even) eigenmodes are, however, compatible with 
having $0 = &2s', on the other hand, from A = 2/, we have 
established that we must have $0 = — $2s> It follows that 
$0 = = 0, i.e. that in the X — 2f mode the insertion 
of spacers that move symmetrically along their bonds leads to 
the immobilization of the black beads. We hasten to note that 
these eigenmodes (involving white and black beads) have the 
same characteristics as the class of eigenmodes discussed in 
Appendix B. 

Given the relative complexity of the above analysis we 
have verified it by considering a series of different, alternat- 
ing black lattices, and by numerically computing the eigen- 
modes and eigenvalues of the corresponding black and white 
systems, obtaining perfect agreement in all cases. 



APPENDIX D 

Here, we show how the results for the eigenmode spectrum 
of generalized dendrimers derived in Refs. II 911 and 12311 can 
be simplified by rewriting them in terms of Chebyshev poly- 
nomials as displayed in Eqs. J60l > and i6l\ of the main text. 
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Following Ref. 1 23], the non-trivial eigenmodes of a gen- 
eralized dendrimer with functionality f c for the central core 
and functionality / for the other junctions belong to two gen- 
eral classes: for the eigenmodes of class (i) the core is mobile 
while for those of class (ii) it is immobile. 

The non-degenerate eigenvalues A r corresponding to class 
(i) can be obtain from Eqs. (6) and (7) of Ref. II23I1 : one has 
namely 



A r = / - 2yJ f - lcos0 r , 
where tj> r is solution of the equation 

sin(g + l)<f> r = - , = 1 sing4> r . 
V J - 1 

Setting x r = cos <\> r these relations are equivalent to 



and to 



\ r = f-2x r y/f=l 



Ug{Xr) = C l Ug-l(X r ), 



which has g solutions x r . In the case 



(Dl) 



(D2) 



(D3) 



(D4) 



(5 + l)/5<l/-/c-l|/V/^l ) 

however, Eq. (ID2> has only g — 1 solutions for real valued 
4> r . Then there appears one more eigenvalue A, which for 
(/ - / c - 1) > is given by Eq. (8) of Ref. Q 



A = /-2V/-lcosh^, 



(D5) 



corresponding to x r = cosh-0 = cos(iip) > 1. Using Eq. 
d!2t with the imaginary angle itp transforms Eq. ( ID4t into 

f — fc — 1 

sinh(<7 + l)tp = - — _L_ smhgip, (D6) 



as given in Eq. (9) of Ref. [23]. In the case that 
{f — fc — 1) is negative, A is given by Eq. (10) of Ref. (2. 



A = f + 2y/f -lcoshtp, 



(D7) 



corresponding to x r = — cosh-0 = cos(tt + itp) < — 1. Using 
Eq. d!2i with the complex angle (n+iip) transforms Eq. iD4\ 
into 

sinh(ff + = - ' ~/ c ~ 1 smhg^, (D8) 

as given in Eq. (1 1) of Ref. l23ll . 

The eigenvalues X r for the modes of class (ii) that are de- 
scribed by Eq. dDlt for real valued angle <\> r obey Eqs. (12) 
and (14) of Ref. |23] with 



sin(m + l)<j) r — yf—1 sin m</> r , with m — 1, g. (D9) 
This can be rewritten with help of Eq.(I12> as 



U m (x r ) = VT- If^m-l(^r) with m = 1, (D10) 
again setting x,. = cos r . The degeneracies are 

A ifm = 5 r 

m ~ I fc{f ~ 2)(/ - l)f- m+1 otherwise. 



This corresponds to Eq. d63l with / c = /. Again, as noted 
in Ref. Q, in the case that m + 1 < m^/J^l, Eq. jD9l 
does not give all the eigenvalues of the modes in terms of real 
valued <fr r . The remaining eigenvalue is given by Eq. iD5i , 
and it corresponds to x r = cosh^ > 1, where now <!/ is 
solution to Eq. (15) of Ref. JH 



sinh(m + — \J f — 1 smb. 771/0. (D12) 



This can again be incorporated into Ea.( ID10> . by making the 
identification x r = cosh?/; = cosz-0. 
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